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Vasily Beskin 



Abstract 

The goal of my lecture is to present the introduction into the 
hydrodynamical version of the Grad-Shafranov equation. Although 
not so well-known as the full MHD one, it allows us to clarify the 
nontrivial structure of the Grad-Shafranov approach as well as to 
discuss the simplest version of the 3 + 1-split language - the most 
convenient one for the description of the ideal flows in the vicinity of 
a rotating black hole. 

1 Introduction 

Axisymmetric stationary flows in the vicinity of a central compact body have 
been studied for a long time in connection with many astrophysical sources. 
Accretion onto ordinary stars and black holes axially symmetrical 

stellar (solar) wind ||, jets from young stellar objects fl , outflow from 
axisymmetric magnetosphere of rotating neutron stars j^,|fj,|^] - all of them 
are flows of the considered type. 

Let me stress that the necessity of taking into account the effects of Gen- 
eral Relativity is not so obvious for many compact sources. For instance, 
one can not exclude that the black hole plays only a passive role in the jet 
formation process, and the effects of General Relativity in this case may 
be inimportant for flow description in the region of jet formation. At the 
same time gravitational effects make, apparently, a noticeable contribution 
to the determination of physical conditions in compact objects. First, this 
is indicated by the hard spectra and the annihilation line observed in galac- 
tic X-ray sources, which are believed to be solar mass black holes. Such 
characteristics are never observed in the X-ray sources which are firmly es- 
tablished to show accretion not onto a black hole but onto a neutron star [|l . 
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Another indication comes from superluminal motion in quasars which may 
be due to the relativistic electron-positron plasma flow ejected along with 
the weakly relativistic jet ||. All this testifies in favour of the existence of 
an additional mechanism for particle creation and acceleration, for which 
the effects of General Relativity may be of principal importance. So, it is 
undoubtedly interesting to consider the flow structure in the most general 
conditions, i.e., in the presence of a rotating black hole. 

There are several reasons why I restrict my consideration to the pure 
hydrodynamic flows. First of all, the hydrodynamical version of the Grad- 
Shafranov equation is not so popular as the full MHD one. On the other 
hand, it contains all the features of the full MHD version in the simplest 
form. In particular, within the hydrodynamic approach one can introduce 
the 3 + 1-split language - the most convenient one for the description of the 
ideal flows in the vicinity of a rotating black hole. 

Thus, in this lecture the basic equations describing a steady axisym- 
metric hydrodynamical flow in the vicinity of Kerr black hole are given. 
Starting with the well-known set of equations describing the nonrelativistic 
ideal flow pxj| , we will go step by step to more complicated cases up to the 
most general one corresponding to the axisymmetric stationary flow in the 
Kerr metric. Finally, several examples will be considered which demonstrate 
how the approach under consideration can be used to obtain the quantita- 
tive description of the real transonic flows in the vicinity of rotating black 
holes. 



2 Ideal Hydrodynamics - Some Fundamental Results 



2.1 Basic Equations 

To start from the very beginning, let us wright down the equations of ideal 
stationary (d/dt = 0) hydrodynamics in a flat space. They are 

• continuity equation 

V-(nv) = 0, (2.1) 



• Euler equation 

VP 

(vV)v = V<^> g , (2.2) 



• isentropy condition 

v • Vs = 0, (2.3) 



• equation of state 

P = P{n,s). (2.4) 
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The last expression can be rewritten in the form 



dP = m p ndw — nTds. 



(2.5) 



Here n (cm -3 ) is the concentration, s is the entropy per particle (undimen- 
sional), w (cm 2 /s 2 ) is the specific enthalpy, m p (g) is the mass of particle, 
T (erg) is the temperature in energetic unit, and c s (cm/s) is the velocity 
of sound. For the polytropic equation of state 



P = k(s)n r , 

used for simplicity in what follows, one can obtain for f ^ 1 



1 

m r . 



dP 
dn 



1 



-Tk(s)n 



r-i 



T 



r- i 

m 



(2.6) 

(2.7) 

(2.8) 
(2.9) 



Now one can make several remarks. 



• The Euler equation (2.2) together with (2__3) and (2__5) can be rewritten 
as the energy equation 

,,2 



V • 



nv 



w 



= 0. 



Now using the continuity equation (|2.]|), one can obtain 

v • VE = 0, 

where 

E = — (c 



v 
~2 



(2.10) 

(2.11) 
(2.12) 



This is the well-known Bernoulli integral. 



• The energy equation (2.10) together with the Euler equation (2.2) can 
be rewritten as a four-component energy-momentum equation 



V„T Q/3 = 0, 



where for ip g = 



^aj3 



nrripV 2 /2 + ne nm p v [v 2 /2 + w) 



P5 lk + nm p v l v k 



(2.13) 



(2.14) 



Here and below the Greek indices a and correspond to four-dimen- 
sional values, while the Latin i, j, and k are three-dimensional. 
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• As a result, hydrodynamics contains five equations over five unknown 
values. 



exercise 

1. Prove expressions (2.1C)-(2.14) 
2.2 An Example - Spherically Symmetrical Flow 

As a most simple but very important example let us consider the spherically 
symmetrical flow. Since the basic (ideal) hydrodynamic equations have the 
form of the conservation laws, one can find 

• the continuity equation 

$ = 2Trr 2 n(r)v(r) = const, (2-15) 

• the entropy equation 

s = const, (2-16) 

• the energy equation 

E= +w(r) + ¥> s (r) = const. (2.17) 

As a result, knowing three parameters $, s, and E, one can determine 
all the ph ysica l characteristics of a flow. Indeed, rewriting the Bernoulli 
equation fl2.17| ) as 

E =^2 r - i + w ^ s )+^{r), (2.18) 

we see that this equation contains only one unknown parameter n. Hence, 
this algebraic equation implicitly determines the concentration n as a func- 
tion of the invariants and radius r: 

n = n(E,s,$;r). (2.19) 

Together with the entropy s, it allows us to determine all the other thermo- 
dynamic functions and the flow velocity v(r ). 

It is necessary to stress that equation ( 2.18 ) contains a singularity on 
the sonic surface. To show this, let us determine the derivative dn/dr: 
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As a result, using the thermodynamic relation one can obtain for the 

logarithmic derivative 

rdn 2v 2 -GM/r 2 — GM/rv 2 N r . „, 

^1 = _ 7T~ = 2 2 = — -i i 212 = ~~n' 2 - 20 



We see that the derivative (2.2C) contains the singularity when the velocity 
is equal to that of sound c* (D — 0). It means that to pass through the 
sonic surface r = r. t the additional critical condition is to be valid: 

CM 

N r {r*) = 2 T =0. (2.21) 

In other words, the transonic flows are determined by two invariants only. 
As shown in Fig. 1, the sonic surface is an X-point on the plane distance - 
velocity. 




Fig. 1. Structure of the spherically symmetric accretion for a given noo and Coo, 
and different stream function $. The transonic flow corresponds to the critical ac- 
cretion rate $ = $ cr (2.34). The curves below the X-point correspond to subsonic 
accretion with $ < $ cr 



exercises 



For spherically symmetrical transonic inflow ( Bondi accretion J|/ ) one 
can determine the Bernoulli integral E through the velocity of sound 
at infinity 



E 



Now using relations (2.1S)-(2.17) and ( 2.2\) , obtain the well-known 
expressions for the velocity of sound c* and the concentration n* on 
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the sonic radius r* 



< = (^3rJ^, (2.22) 
2 \ 1 /( r_1 ) 

n * = 1 5 - 3r ) ( - 2 ' 23 ^ 



5-3r\ GM 



(2.24) 



2. Show that 



, , -4± vio-6r 

^i(n) = ^ , (2.25) 



the sign plus corresponding to accretion and minus - to ejection. 
3. Find that for the spherically symmetric accretion 

• for r 3> r* (subsonic regime) the flow is approximately incom- 
pressible: 

n(r) ~ const, (2.26) 

v(r) cx r~ 2 . (2.27) 

• for r <C r* (supersonic regime) the particle motion corresponds 
to a free fall: 

n{r) cx r~ 3/2 , (2.28) 

v{r) « (2GM/r) 1/2 . (2.29) 

4- Find that for the spherically symmetric transonic outflow (Parker ejec- 
tion fify): 

• Physical parameters on the sound surface r — r*, where again 

r.-f£ (2.30) 
have the following relations to the ones on the star surface r = R 

* - (^H + (^)(»!-T><-> 

/ c 2n V(r-i) 

n* = n R\jprJ ■ (2-32) 
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The (radial) velocity on the star surface is to be 

„2 \ 1/(T-1 



VR 



-R 



(£) ' (2 - 33) 



Although fairly simple, the radial ID flow allows us to formulate here 
several important properties; some of them, as we shall see, retains the 
common properties of the Grad-Shafranov approach. 

• It is possible to pass through the sonic surface in the presence of 
gravity only. Indeed, the nominator N r in ( 2.2d| ) can vanish for D = 



in the presence of the gravity term GM/rv 2 only. 



Solutions (|2.22D -( |2.24j) and fl2.3l|) have singularity for T = 5/3. It 



means that for T = 5/3 the increase/decrease of the velocity of sound 
as a result of adiabatic heating/cooling equals the change of particle 
velocity. As a result, in the nonrelativistic case for T > 5/3 the 
transonic flow is not realized. 

• The transonic problem is two-parametric. It means that to determine 
the transonic flow it is necessary to specify two boundary conditions, 
say, the density m p noo and the velocity of sound Coo at infinity. In 
particular, the accretion rate is fixed: 

/ 2 \(B-3T)/a(r+i) (GM)2 
2$ cr = 4nrlc*n* = 4tt ( — —j g n °°- ( 2 - 34 ) 

On the other hand, for a given and Coo there is an infinite number 
of subsonic solutions with $ < $ cr (see Fig. 1). 

• For a given flow structure the number of integrals is enough to deter- 
mine all the characteristics of a flow from algebraic equations. 

The last property is actually the key point of the approach. Indeed, the 
algebraic relations ( |2.15 )— ( 2.17 ) together with the equation of state allow 



us to determine all the physical parameters of the flow (the flow velocity 
v(r), the temperature T(r), etc.) through the invariants E and s and the 
stream function This property remains true for an arbitrary 2D flow 
structure. But in the general case the flow structure (i.e., the stream func- 
tion <!>(r, 9)) itself is not known. To determine the stream function, the 
extra two hydrodynamic equations are to be used. 



2. 3 Fla t Poten tia I Flow 

To start, let us consider the simplest (and the well-known) case of the flat 
potential flow without gravity. Then one can introduce potential <fi(x, y) by 
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the relation 

v = V0(z,y). (2.35) 
In addition, let us consider the case 

E = const, s = const. (2.36) 

Then the continuity equation V • (nv) = can be rewritten as 

, Vn ■ V(/> , _ s 

V 2 0+ - = 0. 2.37 

n 

Finally, using the Euler equation to determine Vn • W<j) 

f v 2 \ ,Vn-Vc& 



we obtain 



Here again 



j>y) 2 4>xx - ^4>x4>y4>xy + (fyxf^yy 

(V0) 2 L> 



= 0. (2.38) 



£> = -! + -§. (2.39) 



This well-known equation can be found in any textbook (e.g., see JPl] ] 
For us it is necessary to stress here the following properties. 



• To determine c 2 ., equation ( 2.38 ) is to be supplemented with the 
Bernoulli equation. For the polytropic equation of state it can be 
solved for the velocity of sound c s : 

cMr-l)^^) 2 (2.40) 



• Together with the Bernoulli equation, equation ( 2.38 ) contains the 
potential <fi(x, y) and the invariant E only (it does not contain the 
entropy s at all, but s is necessary to determine n) . 

• For n = const {c 2 s — » oo) the equation becomes linear. 

• It is nonlinear in the general case, but linear for second derivatives. 

• The equation is elliptical for a subsonic flow D > 0. 

• The equation is hyperbolic for a supersonic flow D < 0. 

• For a given flow structure (for a given <j), E, and s) all the physical 
parameters are determined by algebraic relations. 
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Equation (2.38) does not contain coordinates x and y. 



The latter property was very widely used in the approach of the hodo- 
graph transformation, i.e., the transformation from a physical plane (x,y) 
to a hodograph plane (v,9), where v x = usin#, v y = vcos8. In this case it 
is possible to int roduc e another potential <f> v (v,6) so that r = W v (j> v . As a 
result, equation ( 2.38| ) can be rewritten as 



d 4> v v d 



= 0. (2.41) 



d6 2 l — v 2 /c 2 dv 2 dv 

This is the linear Chaplygin equation (1902). 

The hodograph transformation approach was the main direction of ex- 
ploration through the XX-th century |l^, |l3). Here I formulate two results 
obtained in this field, to be used in what follows. 

• For the transonic flow it is impossible to solve the direct problem (i.e., 
to determine the flow structure from the given shape of the boundary, 
say, knowing the shape of nozzle or wing) . 

• On the other hand, one can solve the reverse problem. This approach 
is based on the fundamental theorem: the transonic flow is analytical 
at the critical point (the point where the sonic surface is orthogonal 
to the flow line, see Fig. 2) 0. 

Let me comment these two statements. The most visible argument clar- 
ifying the absence of the regular (not iterative) procedure for the transonic 
flow is as follows. As is known, the number of boundary conditions b for 
an arbitrary (not only for a pure hydrodynamical) flow can be determined 
from the following condition JTJ,^5) 

b = 2 + i-cr. (2.42) 

Here i is the number of invariants and a is the number of singular surfaces. 
In pure hydrodynamics the only singular surface is the sonic one. Hence, 
for the transonic flow a = 1. Furthermore, for planar geometry we have 
two invariants, E and s, so that i = 2. Thus, to determine the structure of 
the transonic flow it is necessary to specify three boundary conditions on a 
surface. These may be two thermodynamic functions and one component 
of the velocity. The second (the last in the planar case) component of 
the velocity is to be determined from the solution. But to solve equation 



( 2.38 ), it is necessary to know the Bernoulli integral E = v /2 + w, i.e., 
both components of the velocity on this surface. Hence, in the general case 
even the equation describing the flow structure cannot be formulated. For 
subsonic and supersonic flows a — (so that 6 = 4) and this difficulty is 
absent. 
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v<c, 




sonic surface 



v>c„ 



Fig. 2. Structure of the "analytical nozzle" in the vicinity of the critical point 
x — y — - the only point where the sonic surface is orthogonal to the flow 
line. As the term oc xy in (2A4) is absent (i.e., v y (Q, y) — 0), the plane x — 
corresponds to minimum cross-section of the stream surfaces. The sonic surface 
v = c* has the standard parabolic form x = — [k(T — l)/2c*]j/ 2 



On the other hand, the structure of the transonic flow can be found di- 
rectly by expansion of the solution in the vicinity of the critical point (where 
we put x = y = 0). Indeed, in addition to the invariants E (which deter- 
mines the velocity of sound on the sonic surface) and s (which is necessary 
for the determination of the concentration n) one can specify the x compo- 
nent of the velocity v x (x,0) along x-axis. For our purpose it is enough to 
know the first two terms in the expansion 

v x (x,0) =c* + kx+ ... (2.43) 

Here c\ = 2E(T — l)/(r + 1). As a result, as one can check directly, the 
first terms in the expansion of the potential 4>{x 1 y) look like 

4>(x, y) =c*x+~ kx 2 + -L fc 2 (r + l)xy 2 + -±^k 3 (T+ 1) V + . . . (2.44) 
2 2c* 2Aci 



Knowing all the coefficients in the expansion (|2.43| ), one can restore the 
potential <f> with any precision. 



Thus, in the general case equation ( 2.38 ) cannot be solved directly. 
Let me stress that this property is common; it takes place for the Grad- 
Shafranov equation as well. On the other hand, the planar approach has 
three extra difficulties: 

• It is difficult to consider the case E ^ const, s ^ const. 

• It is impossible to consider the nonpotential flow with Vxv^O. 



• It is impossible to include gravity (which is not planar). 
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3 Axisymmetric Stationary Flow - Nonrelativistic Case 

3.1 Basic Equations 

3.1.1 Stream Function 

Now let us see how a similar procedure can be realized for the axisymmetric 
stationary flows. It means that we assume all the values to depend on two 
variables - r and 9. In this sense the flow remains two-dimensional. But 
now none of the three components of the velocity are equal to zero. Thus, 
axisymmetric stationary flows are more rich than planar ones. 

In the axisymmetric stationary case one can introduce potential <&(r, 9) 
which is connected with the poloidal velocity v p as 

nv p = : — (3.1) 

p 27rrsin<9 v ; 

Such a definition results in the following properties: 

• The continuity equation is valid automatically: V • (nv) = 0. 

• On can check that d<& = nvdS, where S is the area. Thus, the poten- 
tial $(r, 9) is the flux through the circuit r, 9, < tp < 2tt. In partic- 
ular, the total flux through the sphere of radius r is $ tot = 3>(r, 7r). 

• As v- V<& = 0, the relations $(r, 9) = const describe the flow surfaces. 

3.1.2 Integrals of motion 

As earlier, the components /3 = t and (3 = v\\ of the energy-momentum 
conservation law V a T al3 = give 

E = E($)=v 2 /2 + w + (p g , (3.2) 
s = s($). (3.3) 

But as we see, it is now much easier to describe the case when integrals are 
different for different flow surfaces. 

New information appears from the /3 = ip component of the energy- 
momentum equation (or, which is the same, from the ip component of the 
Euler equation) 

V v (^) [v x (V x y)] v + ^ + V^ g = 0. (3.4) 

As we consider the axisymmetric case, all the gradients vanish. The last 
term [v x (V x v)] v can be rewritten in the form of the conservation law 

v V(rsin0t> v ) = 0. (3.5) 
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Hence, in the axisymmctric case the z-componcnt of the angular momentum 

L($) = v v rsin6 (3.6) 
is the third integral of motion. 

3.1.3 Mathematical Interlude - The Covariant Description 

As we are going to generalize our equations to General Relativity, it is 
convenient to rewrite these relations right now in the covariant form. For 
this reason, recall that the flat 3D metric gik (ds 2 = gikdx l dx k ) for spherical 
coordinates x 1 — r, x 2 — 9, and x 3 = ip has a form 

g rr = 1, gee = r 2 , g vv = r 2 sin 2 9, (3.7) 



all the other components being zero. Using now expression (2.14) 



Z? = PS* + [nm p )v k v,, (3.8) 
one can obtain from the tp component of the energy-momentum equation 

1 BP 

\7 k T k = V k (5 k P) + VfcCnmp^) = - — 



r dip 

d_ 

dx k 



+ -^ Zk -(nm p v k v v> ) + T k k (nm p )v t v v - r^nmpje't* = 0. (3.9) 



Here Tj k are the Christoffel coefficients. 

As one can check, the last term in (3.9) vanishes: ^(rimp)?)'^ = 0. 
The first term vanishes because of axisymmetry of the problem (no ip - 
dependence). Using now the continuity equation 



V fe (nm p i; fe ) = — ^Zk(Vgnm p v k ) = (nm p v k ) + T^ k nm p v i = 0, 



^ gd ^^9nm p v) = — k 

(3.10) 

where g = det = grrgeeg^ip , we see that ( |3.9| ) can again be written down 
in the form of the conservation law k T k = (nm p )v • \7v v . Hence, the third 
invariant looks like 

L(4)=V (3.11) 

exercises 



1. Check Eqns. §Jy-§J$. 

2. Is there a contradiction between (3.6) and ( 3.11 ) ? 
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To understand the difference between (3J3) and ( 3.11 ), it is necessary 
to return to the main definitions of the covariant approach. Up to now 
we have dealt with the physical components only. Below in the relativistic 
expressions in Sect. 4 we will mark the physical components by hats, so that 
va — v v is the physical component of the toroidal velocity, and its dimension 



is cm/s. But in the covariant expressions ( 3.9 )— ( 3.11) we encounter other 
objects - contravariant components v l and covariant components Vk- As the 
definition of the vector length (which is square of the physical component) 
has the form v 2 
metric (t 



9ikV "v n 



i 1 ViVk, one can write down for the diagonal 



and the same for the other components. Thus, 



(3.12) 

(3.13) 
(3.14) 

In particular, it means that the dimension of the covariant and contravariant 
components may differ from the dimension of the physical one. Comparing 
now (§T4|) with Q and ( jOl] ), we can understand the difference: (|3.6| ) 
actually involves the physical component of the toroidal velocity while (3.11) 
includes the covariant one. 




3.2 The Stream Equation 

3.2.1 Grad-Shafranov Equation 

To formulate the stream equation (i.e., the equation describing the stream 
function $(r, 8)) it is necessary to return to the poloidal component of the 
Eulcr equation. One can check that together with the definitions of the 
invariants E(&), £(<&), and s(<3?) this vector equation can be written as the 
product of the scalar equation by the vector V<&: [Euler] p = [GS] • V<&. 
For this reason, in many MHD papers the stream equation [GS] = was 
obtained as the projection of the poloidal equation onto the gradient V4>. 
The equivalent hydrodynamic expression has the form 



1 



(V$) : 



(vV)v 



VP 



V^ g 



0. 



(3.15) 



Using now the definitions ( |3.l|) , (|3.2|), and (p.q), we have 



dL 



1. 



-V fc n- V fc $ 



,dE 



T di 



-Att^L— + A^w*n A - 47rWn z — = 



d$ 



d$ 



m p d<& 



(3.16) 



16 



The title will be set by the publisher. 



Here and to the very end 



(3.17) 



so that for the flat metric w = rsin#. 

To clo se the system, i.e., to determine the product (Vn- V<&), the str eam 
equation ( \i.lt ) is to be supplemented with the Bernoulli equation (3^). It 
can now be rewritten in the form (cf. ( [2.18 )) 



E 



(V$) 2 



1 L 2 



, , , , - — +w(n,s) + ip g . 



8TT 2 w 2 n 2 



(3.18) 



As previously, the Bernoulli equation (3.18) contains, besides n, the invari- 
ants E, L, and s, and the stream function $ only. Hence, it again gives the 
implicit expression for n: 



n = n(V®;E,L,s;r,i 



(3.19) 



On the other hand, the implicit Bernoulli equation can be presented in the 
differential form 



where now 



and 



D = -l + -§, 
1 p 



(3.20) 
(3.21) 



V*$ • ViVfc$ , 1 VfeTU 2 , 2 2 2 VfeVg 

-I — • - 47T w n 



-4ttVL 



(V$) 2 
dL V fe $ 
d$ (V$) 2 



(V$) 5 



2n 2 n 2 L 2 



+A-K 2 m 2 n 2 



dE V fc $ 
d$ (V$) 2 



^222 



Vfc^ 2 
zu 2 (V$) 2 
T 1 



dP 

m n n \ ds 



(3.22) 



ds V fe $ 
d$ (V$) 2 



As a result, the stream equation can be written down as UA 



-m 2 X7 k 



(V$) 2 £> 



2L>tu 2 



2 2 2 %' W „ 2 2 g + 1 r 

-An w n — r — rr; 47r n — — — L- — 



+2ir 2 n 



Z?(V$) 2 
2 2 V(w 2 )-V$ 



/i 2 2 2 
-47T TO n 



"D + l T 



L 2 + 47r 2 w 2 n 2 



D d$ 
D + ldE 



(3.23) 



D d$ 



+ 



D m p Dm p n \ ds 



dP 



ds 
d$ 



0. 
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or, in a compact form (cf. p6|), as 



_ ro 2 v ( 1 V fc $ \ _ 4 7T 2 nL^ (3.24) 
/ d$ 

. 2 2 d.E , 2 2 T ds 

+A_-K z vj z n—— - Att zu n — = 0. 

d$ m p d$ 



At first glance, the stream equation (3.23) is much more complicated than 



the planar one (2.38). Nevertheless, one can easily see that these equations 
have many similarities. As ( 2.38| ), the stream equation ( 3.23| ) starts with 



the linear elliptic term and the nonlinear term with an analogous form. The 



third term does not, of course, exist in ( 2.38 ) - it results from non-cartesian 
geometry. But all the other terms are not complification. They allow us to 
include into consideration not only gravity, but a much wider class of flows 
with different invariants onto different flow surfaces. 

In other respects the stream equation is quite similar to the planar equa- 



tion (2.38): 



The stream equation ( p.23| ) is to be supplemented with the Bernoulli 
equation. 



• Together with the Bernoulli equation, equation (3.23) contains the 
potential <&(r, 9) and the invariants E, L, and s only (i.e., it has the 
Grad-Shafranov form). 

• For n — const (c 2 s — > oo), E = const, s = const, and L = the 
equation becomes linear. 

• It is nonlinear in the general case, but linear for second derivatives. 

• The equation is elliptical for a subsonic flow D > 0. 

• The equation is hyperbolic for a supersonic flow D < 0. 

• For a given flow structure (i.e., for a given $) and for the invariants 
E(&), L(&), and s(<i>) all the physical parameters are determined by 
algebraic relations. 

The following point should be stressed. The expression for the denom- 
inator D = — 1 + Cg/w 2 ( 3.21 ) involves the poloidal rather than the total 



velocity. It means that the sonic surface exists when the poloidal, not the 
total velocity becomes equal to that of sound. This fact results from our 
basic assumption: as we consider axisymmetric flows only, the disturbances 
(waves) are to have the same symmetry as well. Hence, the disturbances 
can propagate in the poloidal direction only. For this reason, a singular- 
ity appears at the moment when the particle velocity coincides with the 
velocity of disturbance. 
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3.2.2 Linear Operator tn 2 V fe (n7- 2 V fe ) 

In what follows we shall use the linear operator 

and shall therefore consider it in more detail. First, examine the angular 
operator 

d ( \ d 



* = ™ 6 M(*ToM ■ (3 ' 26) 



It has the following eigenfunctions 



Q = l-cos6>, (3.27) 

Qi=sin 2 6>, (3.28) 

Q 2 = sin 2 6 cos 9, (3.29) 

Q m = (2 l m) , 1 sin 2 (cos*), (3.30) 



and eigennumber values 



q m = -m(m + l). (3.31) 



Here P m are Legendre polynomials and prime means their derivative. As a 
result, for the full operator C we have the following set of eigenfunctions 



1. m 



$^ = r 2 sin 2 6 - homogeneous flow, 



• $^ = sin 2 9/r dipole flow. 

2. m = 2 

• $^ — r3 sm2 ^ cos ^ ~ zer0 point, 

• $^ = sm2 0cos8/r 2 - quadrupole flow. 

3. ... 

At first glance, this point is absolutely clear so that it is impossible 
to encounter here any trouble. Nevertheless, it is not so. Indeed, let us 
consider the eigenfunctions corresponding to m — 0. The first one is clear: 
it is the function 

§ { q> = 1-cos6 (3.32) 

which describes the spherically symmetric accretion/ejection. By the way, 
this harmonics alone determines the accretion/ejection rate as for all the 
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other eigenfunctions with to > we have $ m (r, ir) = 0. The uncertainty is 
due to the second eigenfunction 

$( 2 > = r (l - cos (9) (3.33) 

for which <J> (r, 7r) ^ const. It means that this harmonics can be realized 
only if in the volume (not in the origin or at infinity) there are the sources 
or the sinks of matter. In all other cases the second eigenfunction for to = 
is to be dismissed. 



3.3 Examples 

3.3.1 Bondi-Hoyle Accretion 

As a first example we consider the accretion onto a moving gravity center 
(Bondi-Hoyle accretion Q). This is a classical problem of modern astro- 
physics. The nature of Active Galactic Nuclei and Quasars, the nature of 
jets, the activity of some galactic X-ray sources are believed to be associated 
with the accretion of a gas onto compact objects - neutron stars and black 
holes [§,17, Nevertheless, only a few exact solutions describing the 



accretion flow (even for the simplest adiabatic case) are now known, e.g., 
the solution for the spherically symmetric flow we have already discussed. 

To construct the nonsphcrical solution, one may assume that a small 
perturbation of a spherically symmetric flow cannot change strongly the 
structure of the accretion pC|| . So, it is possible to seek the solution of the 
stream equation as a perturbation of the spherically symmetric solution. 
First of all, let me recall the main results of the qualitative theory. It is more 
convenient to perform the calculations in the reference frame moving with 
the gravity center with a velocity . Comparing the Bondi accretion rate 
2$ = AnrlrivC* ~ (GM ) 2 n 00 /c^ (^.22[)-(|2~2|) with the flow $ - irR^n^Voo 



captured within the capture radius R c , one can evaluate the R c value as 

R c ~e- 1/2 r* : (3.34) 

where 

ex = — . (3.35) 

Coo 

Hence, for e\ 1 the capture radius is much larger than the sonic one, 
and we can assume that for r <C R c the flow structure is similar to the 
spherically symmetric accretion. Thus, one can seek the solution of the 
stream equation ( 3.23| ) in the form 



$(r,0) = $ [1 - cos 6 + £l f{r, 6)}. (3.36) 
For a nonmoving gravity center we return to the spherically symmetric flow. 



As the stream equation (3.23) now contains i = 3 invariants, so that 



b = 2 + 3 — 1 = 4, it is necessary to specify four boundary conditions, say 
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1- Voo — const, 

2. v v — (and hence L = 0), 

3. Sqo = const, 

4. £ 00 = (4/(r + i). 

In the last relation we neglect the terms ~ e\. 

As a result, the stream equation can be linearized: 



n d2 f £l m^-n • a d ( 1 GM\df 



96» Vsin6» 90 y ' ^ V?- c s 2 r 2 y cV °" *' 3 ' 37 ' ) 



This equation has the following properties. 
• It is linear. 



The angular operator coincides with Cg (3.26) 



• As all the terms contain a small E\ value, the functions D, c s , etc. can 
be taken from the zero solution. 

• As for the spherically symmetric flow the functions D, c s , etc. do 



not depend on 6, the solution of equation (3.37) can be expanded in 
eigenfunctions of the operator Cg. 

Thus, the solution can be presented in the form 

oo 

f(r,d) = 9m{r)Qm(0), (3.38) 

m=0 

the equations for the radial functions g m {r) being 



r 2 D f^L + ( 2r - ^f) % + m(m + 1)(D+ l)g m = 0. (3.39) 
or* \ Cg / dr 

As to the boundary conditions, they can be formulated as follows: 

1. No singularity on the sonic surface (where r 2 N r = 2r — GM/c 2 = 0, 
D = 0) 

<?m(r*)=0. (3.40) 

2. A homogeneous flow $ = irrioaVoar 2 sin 2 9 at infinity, which gives 
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As a result, the complete solution can be presented in the form 
$(r, 9) = $q [1 - cos 9 + em (r) sin 2 9)} , 



(3.42) 



where the radi al function gi(r) is the solution of the or dinar y diffe renti al 
equation ( |3.39 ) for m = 1 with the boundary conditions ( 3.4C ) and ( |3.41 ). 

At the present level of personal computers it means that we succeed in 
constructing the analytical solution of the problem in question. It allows us 
to obtain all the information concerning the flow structure. In particular, 
the sonic surface now has non-spherical form: 



r„(0) =r. l + 2ei 



1 



D 



&2 cos 9 



(3.43) 



where the numerical coefficient &2 = t*9\{i'*) can be obtained from the 
solution (for more details see p0|]). As one can see from Fig. 3, the analytical 
solution is in full agreement with the numerical calculations |21| although 
the parameter e\ = 0.6 is not too small. Finally, the condition 



.9o = 



(3.44) 



(to obtain which an additional consideration is necessary) shows that the 
accretion rate is not changed in this approximation. 



05 (M/cl) 
0.2 



0.1 




z (MK) 



Fig. 3. Flow structure and shape of the sonic surface for T = 4/3, si =0.6 @. 
Labels on the curves denote the values of <E>/$o, and dashed curves indicate stream 
lines and the sonic surface obtained numerically in pl[ 



In connection with this solution, it is necessary to clarify one point. 
As one can easily see, outside the capture radius, our main assumption - 
smallness of the disturbance of the spherically symmetric flow - is not valid. 
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Nevertheless, the constructed solution is correct. This beautiful property is 
connected with the already mentioned fact that for an (approximately) con- 
stant concentration n the stream equation becomes linear. But as we know 



from the Bondi accretion (see (2.2q)), far from the sonic surface r ^> r„ the 



flow density is actually constant. The same is true for the homogeneous 
flow. As a result, for R c 3> r*, i.e., for the case £1 <C 1 under consideration, 
in the vicinity of and outside the capture radius (where the "disturbance" 
~ £igi(r) becomes of the same order as the zero approximation ~ 1), the 
stream equation is linear. As a result, the sum of two solutions, homoge- 
neous and spherically symmetric, remains a solution as well. 

3.3.2 Ejection from a Slowly Rotating Star 

Another interesting nonrclativistic example is the transonic ejection from a 
slowly rotating star p2],p3|. It is necessary to stress from the very beginning 
that this example is only illustrative because in reality an important role is 
played by the radiation pressure which cannot be included into consideration 
within our approach. Nevertheless, the analysis of this problem helped us 
to clarify some important features of the Grad-Shafranov approach p4| . 
As a zero approximation we consider the well-known Parker transonic 



outflow (2.30)— (2.32). It means that we assume all the parameters of the 
spherically symmetric outflow (the sonic radius r* , the velocity of sound on 
the sonic surface c*, the radial velocity vr on the star surface r = R, etc.) 
to be known. As before, we will seek the solution in the form 

$(r,0) = $o[l-cos0 + 4f(r,0)], (3.45) 

where the small parameter is now 



GM 



(3.46) 



Here is the angular velocity of a star. 

The problem under consideration needs all the i = 3 invariants. Hence, 
6 = 2 + 3 — 1 = 4, and it is necessary to specify four boundary conditions 
on the star surface r = t\r(0) which now differs from the sphere 

r R (9) = R[l + e 2 2 p(6)}. (3.47) 

Here we introduce the dimensionless parameter p(9) ss 1. 

Let me stress that at first glance there is a disagreement as we add one 
degree of freedom (the toroidal rotation with v v ^ 0) while it needs two 
extra functions in comparison with the spherically symmetric outflow. This 
question will be clarify below. 

It is important that for small £2 four boundary conditions can be deter- 
mined through real physical parameters on the star surface, e.g., through 



Bcskin: 2D Transonic Hydrodynamics 



23 



two thermodynamic functions (say T and n) and two components of the 
velocity (say v r and v v ). As a result, one can express these boundary 
conditions through four dimensionless functions t(6), f](9), u)(9), and h{9): 

T(r R ,9) = T R [l + e 2 r(9)}, (3.48) 
n(r R ,9) = n R [l+4ri(0)], (3.49) 



GM\ 



1/2 



v v (r R ,6) = £ iy—) u(0)sm9, (3.50) 

v r (r R ,9) = v R [l + e 2 2 h{9)}. (3.51) 
Here ui(9), determined as 

Q(r R ,9) = ih){6), (3.52) 

describes the differential rotation of the star surface. 
Using now the thermodynamic relation 

1 dT dn 

d ^T-I---' ^ 



one can obtain for the first invariant s(9) 

1 



Ss(9) 



r - 1 



t(9) - r,(9) 



(3.54) 



Accordingly, two other invariants can be determined through the boundary 
conditions as well 

CM f 2 F T 

5E{9) - elv 2 R h{9)+el w 2 (9) sin 2 6 + -2—— t (6) + 6^.55) 

Lti 1 — 1 nip 

L 2 (9) - £ 2 R 2< ^-w 2 {9) sin 2 9. (3.56) 

Here we use the following expression for the disturbance of the gravitational 
potential on the star surface 

6<p e (r R ,9)=s 2 ™p(9). (3.57) 

It is of great importance that the possibility of taking the next step, i.e., 
writing down the Grad-Shafranov equation, is connected with the simple 
geometry of the zero approximation. Since in the zero approximation the 
stream function is <f> = 1 — cos 9, i.e., it is the function of 9 only (and as all 
the derivatives di?/d<!>, dL/d<I>, and ds/d$ as well as L itself vanish for a 
nonrotating outflow), one can use the relation d<i> = $o sin0d#. It allows us 
to determine the derivatives d£?/d$, dL/d$, and ds/d$ in the whole space 
and not only on the star surface. 
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As a result, the stream equation can be written as 



-e 2 <f> D^ - 4$oCD + l)sin0- 
or z r z 



8 



1 df 



80 \ sin 86 



df 
Or 



-A-n 2 n 2 r 2 sin 2 6(D + 1)— + 47r 2 n 2 (L> + 1)L^ 
d<I> d<f> 



2 2 cos " r 2 i a 2 2„2 • 2 /i 

-An n ^ — L + Air n r sin 

$n sin 



r r- 1 



(D + l)— + 

m n r 



(3.58) 
ds 



Here N r = 2/r - 47r 2 n 2 r 2 GM/$ 2 . 

The properties of this equation are the same as before. 

• It is linear. 



• The angular operator coincides with Ce ( 3.26 ). 

• As all the terms contain a small s 2 value, the functions D, c s , n, etc. 
can be taken from the zero approximation. 

• As for the spherically symmetric flow the functions D, c s , n, etc. do 
not depend on 9, the solution of equation ( 3.58| ) can be expanded in 
eigenfunctions of the operator Cg. 



Hence, we can again seek the solution in the form 

oo 

f(r,6) = J2 9m(r)Qm(0). 



(3.59) 



m— 



Introducing the dimensionless variables 



n c 
u = — , a = - 

C: 



(3.60) 



one can write down the following ordinary differential equations describing 
the radial functions g m (r) 



[ 1 — x au ) — h 2{ x u )— h m(m + l)a; au g n 

1 x ' ' 



da; 2 



da; 



4 4 

x au 



,■2 



Am. — ^-U 2 



6 4 

u m x au 



--v m x 6 a 2 u A 



r- 1 



(3.61) 



v m x 2 au 2 , 
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where n m , A m , a m , and v m values arc defined as the expansion coefficients 
AE 00 ^ 

m— 

00 

= 4c>t J2 X rnQ m {0), 



d9 



cos 9 
sin 2 



-L 



m=0 



L AL 

^e~AB 

ds 



4 C l r l K rnQm{0), 
m=0 

00 

£2 ^ v rn Q m (&). 



(3.62) 
(3.63) 
(3.64) 
(3.65) 



Finally, the functions a{x) and corresponding to the spherically 
symmetric flow for the polytropic equation of state ( |2.6| ) are related as 
a = u T ~ x and the function u(x) due to ( |3.2C )-( 3.22| ) can be found from the 
ordinary differential equation 



dit 1 — x 3 u 2 

dx x 1 — x 4 au 2 



with the boundary conditions (cf. (2.25)) 



»-='- l's 



4 + vio-er 
f+1 ' 



(3.66) 



(3.67) 



As to the boundary condition to the set of equations ( |3.61 ), they are 
quite similar to the Bondi-Hoyle accretion (for more details see [p4[): 



• Condition on the star surface. As 



d$ = 2nr 2 nv r sin 6A6 = 2TrR 2 n R v R [l + e2(r) + h + 2p)]sm6A0, (3. 
we have 

(2m)! 



g m (R/r* 



2 m {m + l)!m! 



+ 2 Pm ). (3.69) 



Here r? TO , /i m , and p m are expansion coefficients in Legendre polyno- 
mials, e.g., r](9) = Y, m VmP m (cose). 

• The absence of singularity on the sonic surface Ng = 0. This condition 
gives 

(2m)! 



£ 2 ff™(l) 



2 m (m + l)!m! 



(<J£7) m . (6L/ sin 2 0) r 
- (os) m - 



c-t 2c^rJ 

(3.70) 

where again (. . .) m means the expansion in Legendre polynomials 
which can be found from expressions ( |3.54 )-( 3.56| ). 
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As a result, equations (3.61) taken together with the boundary conditions 



(3.69) and ( 3.70| ) allow solution of the direct problem, i.e., determination of 



the flow structure from the physical boundary condition on the star surface. 
Here it is necessary to stress two important points: 

1. The reason we succeeded in writing the regularity condition Ng = 
on the sonic surface (and hence to solve the direct problem) is again 
a very simple geometry of the zero approximation. In particular, in 
the problem in question the position of the sonic surface can be taken 
from the spherically symmetric solution. In the general case it is not 
so and it is impossible to write down the condition Ng — using the 
known functions (SE) m , (8s) m , etc. on the star surface. The position 
of the sonic surface is unknown and is to be found from the solution. 

2. We can now understand the nature of the "extra" boundary condition. 
The point is that for m — one can take, as was already stressed, one 
fundamental solution only, namely, go = const. Another fundamental 
solution is unphysical. Hence, there is an additional relation 

g (R) = 5b(r„). (3.71) 

This relation determines ho which, as we see, is not a free parameter. 
In other words, we are not fully free in the boundary condition h(9): 



its zero harmonics is to be found from relation (3.71). But as was 
demonstrated earlier, it is go that determines the ejection rate. Hence, 
the ejection rate is a function of three parameters only, namely, zero 
harmonics of two thermodinamic functions rjo and To, and v v . For 
a spherically symmetric flow v v = 0, and we return to two functions 
which determine the ejection rate. As to higher harmonics with m > 0, 
they are free, and to determine them it is necessary to know four 
functions on the star surface. Thus, the spherically symmetric case is 
degenerate and it is necessary to be very careful when extending its 
properties on the 2D flows. 

At the end of this section let me formulate some results which can be 
obtained under the following simplified assumptions: 

• Almost the whole of the star mass is in its center, i.e., ip s = —GM/r. 

• No differential rotation, i.e., u>(6) = 1. 

• Von Zeipel law: T(R, 9) oc gta', where <p c g = ip g + L 2 /r 2 . 

• No meridional convection, i.e., Vg(rR, 9) = (it means that we specify 
here vg(rn,9) instead of v r (rn,9) , i.e., the coefficients ho, hi, etc. are 
to be found from the solution). 
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exercises 



1. Find that the distur bance s of the star radius p(9) in (S^Jj) and the 
temperature t(9) in (3.46) have the form 



P{0)= 2 sin 2 



7-09) 



■ sin 



(3.72) 



2. Show that the only nonzero terms in the expansion 1 3.5 i ) correspond 
to m = and m = 2. the expansion coefficients in (3.6i )-( 3.65\ ) being 



<T 2 

A 2 



2^A^ + i4-3% 2 , 



r 2(r-i) 



2 c 2 



r r r 

2 — , k 2 = 4— , v 2 = -- — -. 

r* r* 1—1 



(3.73) 
(3.74) 



and (To, vq = 0. Remember that h(9) = ho + h\ cos + /i2-P2(cos ( 



As a result, solving the stream equation (3.61) for to = 2, one can find 



that 



1. The ejection rate can be presented as 

q 2 r 3 



$tot = 2$n 



1 



GAI 



(i + V) 



Here /io can be obtained from relation ( 3.71 ) (see Table 1) 
, 1 2n/R-R/u 



G 



3 l-vl/cl 



(3.75) 



(3.76) 



As we see, the rotation increases the ejection rate. 
Table 1 



model 


l + h 


hi 


92 


bo 


62 


r*/R= 1.1, T = 4/3 


2.9 


-0.8 


-0.40 


2.2 


-0.41 


r*/#= 2.0, T = 4/3 


3.2 


-3.3 


-0.71 


2.5 


-0.47 


r„/i? = 10, T = 4/3 


8.0 


-56.0 


-2.18 


5.7 


-0.92 


r*/i?= 1.1, r = 1.1 


1.7 


-0.8 


-0.15 


1.8 


-0.37 


r*/R= 2.0, T = 1.1 


2.1 


-2.3 


-0.18 


2.2 


-0.40 


r*/H= 10, r = 1.1 


7.4 


-26.0 


-0.40 


7.2 


-0.58 
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2. Far from the sonic surface r»r, the stream function has the form 
*(r, 

— = 11 — cos c i i — 

v ' GM 

n 2 R 3 



lim 

r — >-oo <J>n 



Q 2 R 3 

(1 - cos0) + --^-(1 + M(l - cosf 



GM 



■q 2 sin 0cos9, (3.77) 



where the coefficient (72 is tabulated in Table 1 as well. As there is no 
r-dependence, the outflow becomes purely radial at large distances. 



3. Accordingly, the asymptotic expression for the concentration n 
the form 



has 



lim 



n(r, 9) 



- b o + 7; ^77-^2 (3 cos 2 6> - 1) 



GM 



2 GM 



where 



(3.78) 

2GM/R. As one can see from Table 1, b 2 < 0. It 
means that the rotation results in the appearance of a dense disk in the 
equatorial plane; this result is well-known |23|, but previously it was 



obtained by numerical calculations only. Negative q 2 values in (3.77) 
demonstrate that for e 2 > 1 the most mass outflow is concentrated in 
the vicinity of the equatorial plane as well. 



4 Axisymmetric Stationary Flow - General Relativity 

4.1 Basic Equations 

We shall show how the Grad-Shafranov approach can be applied to the 
axisymmetric stationary flows in the vicinity of a rotating black hole. Re- 
member that the main difficulty of the General Relativity is the necessity to 
work with four-dimensional objects. As a result, we cannot use our three- 
dimensional intuition in considering the relativistic processes. 

But there is a convenient language - 3 + 1-split - which allows work 
with three-dimensional vectors even in General Relativity [ p5| . One can 
find the detailed introduction into this approach in the book "Black Holes. 
The Membrane Paradigm" by K. Thorne, R. Price, and D. Macdonald [ p7| . 
The main idea is that for stationary metric the proper time r and "the 
time at infinity" t are in one-to-one correspondence. It allows the time t 
to be separated from spatial coordinates x l (i — 1,2,3). As a result, all 
the equations can be written in the simple 3D form, their physical meaning 
remaining clear. Below I shall give the main relations of this approach 
(see [^6| as well). 
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4.1.1 Kerr Metric 

The Kerr metric is the metric of a rotating black hole. In the Boyer- 
Lindquist coordinates t, r, 9, and tp it has the form 

ds 2 = -a 2 dt 2 + g ik {dx l + ^dt)(dx k + (3 k dt), (4.1) 

where 

a=|VA (4.2) 
(a is the lapse function or the red shift), and 

0^ = ^=0, i u> (4.3) 

(oj is the Lense-Thirring angular velocity; remember that (3 V is the con- 
travariant component). Finally, it is important that the 3D metric gik in 
(hi) be diagonal 

P 2 

9rr = -^, gee = p 2 , g vv = w 2 . (4.4) 

Here M and a are respectively the black hole mass and the angular momen- 
tum per unit mass (a = J/M), and we introduce the standard notation 

A = r 2 +a 2 - 2Mr, p 2 = r 2 + a 2 cos 2 6, 

X 2 = (r 2 + a 2 ) 2 — a 2 A sin 2 9, m = -sm9. (4.5) 

P 

Units where c = G = 1 are used below. 

This metric has the following properties. 

• The Kerr metric is axisymmetric and stationary. It is precisely what 
is needed for using the Grad-Shafranov approach. 

• The Kerr metric is two-parametric, i.e., it depends on two parameters: 
the mass M and the rotation parameter a. 

• It transforms to the Schwarzschild metric for a nonrotating black hole: 
g rr — or 2 , g gg — r 2 , g vv — r 2 sin 2 9. Here a 2 — 1 — 2M/r. 

• In the limit r ^S> 2M Boyer-Lindquist coordinates coincide with the 
spherical ones: g rr = 1, ggg = r 2 , g vv — r 2 sin 2 9. 

• Lapse function a 

— The lapse function a describes the lapse (the red shift) between 
the proper time r and the time at infinity t: dr = add. 
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The condition a = determines the position of the horizon 



r g = M+ \JM 2 - a 2 . (4.6) 

— Boyer-Lindquist coordinates do not describe the space-time in- 
side the horizon; for r — r g the metric has a coordinate singular- 
ity. 

• Lense-Thirring angular velocity u 

— The Lense-Thirring angular velocity uj corresponds to the proper 
motion of the space around a black hole. 

— By definition, Sljj = w(r g ) is the black hole angular velocity (does 
not depend on 9). 

— uj oc a, r^H^g = a/(2M). 

• Convenient reference frame - ZAMO 

— ZAMO (Zero Angular Momentum Observers) [^?J are located at 
a constant radius r = const, 9 — const, but they rotate with the 
Lense-Thirring angular velocity dip/dt = uj. 

— For ZAMO the four-dimensional m etric g a p is diagonal, the spa- 
tial 3D metric coinciding with (4.4). 

— No gyroscope rotation in the local experiment. 

To clarify the physical meaning of the values a and to, let us consider 
the motion of particles in the gravitational field of a rotating black hole. 
Then, one can rewrite the four-dimensional equation of motion 

d 2 x a n dx$ dx~f , , 

^ + r ^i7" ^ 

in the simple 3D form 

d Pi _ m p m p v k 

d7 ~ vT^? TT^ 1 ( 8) 

where 

g = — Va, (4.9) 
a 

H lk = -Vifa- (4.10) 

a 

Remember that 
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• The Greek indices a, 0, and 7 are four-dimensional, while the Latin 
i, j, and k are three-dimensional. 

• t is the proper time, and all three-dimensional vectors are measured 
by ZAMO. 

• V,- means the covariant derivative in the three-dimensional metric 



(4.4) 



In a weak gravitational field, i.e., far from a black hole there is very nice 
analogy between the gr avit ational and electrodynamic equations. Indeed, 



the equation of motion (4.8) can be rewritten as 



where 

g=-Va, H = Vx/3, (4.12) 

a and /3 playing the role of the scalar and vector potentials, respectively. 
Moreover, the Einstein equations in a weak gravitational field are quite 
similar to the Maxwell equations 

V g = -47Tftn, (4.13) 

V x g = 0, (4.14) 

V H = 0, (4.15) 
V x H = -167rp m v. (4.16) 

In other words, the gravitational field g is analogous to the electric field 
while the new (so-called gravitomagnetic) field H - to the magnetic one 
which is proportional to the angular velocity of a black hole. The sources of 
the gravitoelectric field g are masses, and the sources of the gravitomagnetic 
field H are mass currents. 

For example, for a rotating sphere with mass M and angular momentum 
J the fields outside the sphere are ]27j 

M 

g = -^e f , (4.17) 
H = 2 J -y^) , (4.18) 

i.e., the rotation induces a dipole gravitomagnetic field around the rotating 
body. The appearance of an additional gravitomagnetic force is the most 
important consequence of the black hole rotation. 

To summarize, the 3 + 1 language allows the description of physical 
processes in a clear 3D form. If we simultaneously use ZAMO as a reference 
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frame, no expressions will contain extra terms. In a sense, ZAMO is an 
inertial frame, at any case in the ip direction. As a result, within the 3 + 1- 
split 

• All three-dimensional vectors are to be determined from the local ex- 
periment by ZAMO. 




• All the calculations are to be made in the 3D diagonal metric (4.4), 
e-g-, 

V-A = ~(V-9A% (4.19) 

y/gee e § V#w e v \ 
d/dO d/dip . (4.20) 

y/geeAg ^/g^pA,p J 

• All the vector relations remain the same as in the flat space, e.g., 
V x (Va) = 0, V • (V x A) = 0. 

4.1.2 Thermodynamics 

On the other hand, all thermodynamic functions within the 3 + 1 approach 
are determined in the comoving reference frame, which is the only invariant 
one. For this reason, one need no think about the transformation connected 
with different reference frames. Actually, there is the only complification: 
in the relativistic case one should work with the relativistic enthalpy /i 
including the particle rest mass 

\i = — m m p c 2 + m p w + . . . (4-21) 

Here p m is the internal energy density. For the polytropic equation of state 
P = k(s)n r (|J) we have (c = 1) 

M = m p + -^—k{ S )n r ~\ (4.22) 

cl = U?f) =-k{s)r?-\ (4.23) 

Finally, the relativistic energy-momentum tensor has a symmetrical form 

p ( £ S\_/ ( Pm + /V) 7 2 (p m + P) 7 U \ 

IS T lk ) I (p m + P) 7 u (p m + P)u*u k + Pg* k ) 



Remember that 7 is the Lorentz- factor of a flow measured by ZAMO. 
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Now, using the relativistic version of the energy- momentum conservation 
law V a T af} = 0, one can obtain (25| 



-(/3V)e = -iv • (a 2 S) + H lk T lk , 
a or 

du> 



a 



a 



<9x fc 



(4.25) 
(4.26) 



Here the additional terms in the energy (4.25) and momentum (4.26) equa- 
tions are due to the gravitomagnetic forces. 



4.1.3 Stream Function, etc. 

As in the flat space, one can introduce the stream function <I>(r, 9) through 
the poloidal component of the four- velocity of the flow u p 



anu r 



p 27rtn 

It means the following relations for the physical components 



1 



anuf 



2ttzu 
1 

277-07 



(V*) f . 



The definition (4.27) gives the continuity equation in the form 

V • (emu) = 0. 



(4.27) 

(4.28) 
(4.29) 

(4.30) 



Let me clarify the extra factor a in (4.30). The point is that the 3D conti- 
nuity equation ( 4.30| ) results from the 4D one 



d 



(Vs^ 3 ), 



(4.31) 



where gu is equal to a 2 . 

As a result, using the definitions (4.24) and (4.27), one can rewrite the 
energy equation ( 4.25| ) and the ip component of the momentum equation 
(gTggh as 



u ■ V(a/i7) + fiu^u ■ Vw = 0, 
u ■ V(//%) = 0. 

Hence, two integrals of motion can now be present as 



E(<f>) 



(4.32) 
(4.33) 



(4.34) 
(4.35) 
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exercise 



Obtain expressions (4-3i)-(^.3L). 



Expre ssion s ( 4.34 ) a nd ( 4.35 ) are the extension of the nonrelativistic 
relations (3^2) and (3J3) to the case of a rotating black hole. Indeed, for 
uj = we have for Bernoulli integral 



1 v" 



E = 7/ia « ( 1 + J (m p c + m pW ) ( 1 - ^, 



GM 



TOpC + TO p [ — + W + 93 g 



As to the third invariant, we again have 

s = s(4>). 



(4.36) 



(4.37) 



4.1.4 Grad-Shafranov equation 

Using now three invariants E, L, and s, one can write the poloidal compo- 
nent of the relativistic Euler equation |31fl 



w 



1 1 2 

— fiwylwusW w H — V a a = 
a a 



(4.38) 



(here the indices a and b are r and only) as [Euler] p = [GS]V$, where the 
stream equation [GS] = now looks like 



Here 



-M 2 



avj 2 V k 



J_ vfe$ \ , V a -$V 6 -<I>V (1 V ( -,<T» 



(V$) 2 L» 



2(V$) 2 L> 



(4.39) 



32tt 4 9 
M 2 d$ 

64tt 4 



k7 2 (£ - loL) 2 - a 2 L 2 ] - 16n 3 a 2 zu 2 nT-^- = 0. 

d$ 



F 



[w 2 (E - ujL) 2 - a 2 L 2 - n7 2 «V] 



and we introduce the thermodynamic function 

M 2 = 



(4.40) 



(4.41) 
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Next, the operator V' fc acts on all the variables except A4 2 , and d/d<& on 
the invariants E(&), L(Q), and s($) only. Finally, now the denominator D 
is 

£> = -L+^-iL. (4.42) 
Here the physical component of the poloidal four-velocity u p can be found 



from (4.27). In a compact form the relativistic stream equation looks like 



32tt 4 d 
M 2 d¥ 



[w 2 {E-wL) 2 -a 2 L 2 ] 



-mr'a^^nT— = 0. 

d$ 



(4.43) 



The hydrodynamical version of the Grad-Shafranov equation in the Kerr 
metric was first formulated in |28| using four-dimensional notation and 
in [^9j within the 3 + 1-split language. 

As before, the Grad-Shafranov equation is to be supplemented with the 
Bernoulli equation (j 2 — 1 + + u 2 ) which can now be written as 

(E - ujL) 2 



2 2 

a (i 



M 4 



64ir 4 vu 2 



(4.44) 



Remember that in ( 4.39| ) and ( 4.44 ) the relativistic enthalpy \x is to be 
considered as a function of M 2 and s: n = fj,(A4 2 ,s). In the general case 
the appropriate differential connection is |2D| 



AM 2 



1 - c? 



dP 

ds 



T 



As. 



(4.45) 



In particular, as in the nonrclativistic case, the Bernoulli equation (4.44) 
determines M 2 in the implicit form through the flux <!> and three integrals 
of motion: M 2 = 7W 2 (V$; E, L, s). 

Finally, it is convenient to use another form of the poloidal four- velocity 



E 2 -a 2 L 2 /w 2 



a 2 ^ 2 



a yr 



(4.46) 



We see that u p — * oo as a 1 at the horizon. As was already stressed, 
it results from our choice of the reference frame which has a coordinate 



singularity for r = r g . Relation ( |4.46 ) suggests a very important conclusion. 



The flow in a close vicinity of a black hole is to be supersonic (u p > c s ). 
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4.2 Examples 

A. 2.1 Exact solutions 

1. Spherically symmetric accretion 

If the flow velocity at infinity is equal to zero (700 = 1) and ther- 
modynamic conditions are homogeneous, then E = = const and 
s = Soo = const. So, the two thermodynamic functions at infinity 
determine two integrals of motion E and s. Finally, for a spherically 
symmetric flow one can put L — 0. Under such conditions the stream 



equation ( 4.39 ) has a trivial solution $ = $ (1 — cos#), and the ac- 
cretion rate 2$o is to be determined from the regularity conditions on 
the sonic surface r = r*. 

As a result, we have the following expression for the sonic radius 

r. -£■('+»), (4.47) 



so that for « 1 we return to the nonrelativistic expression ( |2.30 ). 
As to the relation between c\ and c 2 ^, it can be found from the con- 
dition n 

1 = j^MZil+M) (448) 

[i-<4/(r-i)] 2 [i- c 2/(r-i)] 2 .(i- c 2)- 

In the limit < 1 we return to the well-known relation ( |2.22| ). Next, 
the values M\ and /x* on the sonic surface are 



/j \ V(r-i) / r _ 1 _ „ 2 x (2-r)/(r-i) 

r - 1 - c 2 

= _ T - ( 4 - 5 °) 

Thus, the accretion rate can be written as 2m p $o, where we now have 
$0 = 8ir 2 rlEc*/Mi 



exercise 

Show that in General Relativity a transonic flow takes place even for 
T = 5/3. 

Finally, in the supersonic region r <C we have 
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In particular, at the horizon 



-3(r-i)/2 

(4.52) 



c s(^)-Y77f-r( c *) 5_3r - ( 4 - 53 ) 

Hence, for c 2 w -C 1 the velocity of sound remains small (c s <C 1) 
up to the horizon. 

It is necessary to stress that an accretion onto black holes differs suffi- 
ciently from the nonrelativistic case. The point is that all the subsonic 
trajectories taking place for an accretion onto ordinary stars (see Fig. 
1) have unphysical singularity v(r — > r g ) = 0, n(r — > r s ) — oo on 
the horizon. In other words, to support the subsonic flow the infinite 
gravitational force in the vicinity of the horizon is to be balanced by 
the infinite pressure gradient. Hence, the only physically reasonable 
regime of the accretion onto a black hole is the transonic one. 

2. Accretion of a dust (P = 0) 

For a dust, the flow lines must coincide with trajectories of particles 
freely moving in the gravitational field of a Kerr black hole. For the 
case of the motion with L = (u v = 0) and zero kinetic energy at 
infinity 7oo = 1, such trajectories are "straight lines" 9 = const for an 
arbitrary rotation parameter a pi] ]. Moreover, for P = the density 
of the flow lines can be arbitrary as well. In other words, the arbitrary 
function 

$ = $(6>) (4.54) 

must be a solution to the stream equation. In particular, it means 
that the accretion rate is free. It is not surprising, for the flow is 
supersonic in the entire space. 



exercise 



Using the Bernoulli equation ( 4-44 ) an d the compact form of the rel- 
ativistic stream equation (£.4&), check that for E = fj, = const, L = 0. 
and s = the arbitrary function is a solution. 



3. Accretion of a gas with c s = 1 



1 we have D 1 = 0. Hence, for 
const, L = 0, and s = const the stream equation becomes linear 

A a 2 $ 



As one can see from (4.42), for c. 
E 



dr 2 



1 • n d 



1 <9$ 



sin 9 89 



0. 



(4.55) 
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As a result, the solution can again be expanded in eigenfunctions of the 
operator Cg. For example, for a moving black hole one can obtain (3TJ] 

$ = $o(l - cos6>) + im^v^r 2 - 2r g r) sin 2 9. (4.56) 

Here the accretion rate 2$o is arbitrary for another reason - the flow 
remains subsonic up to the horizon. 



4.2.2 Bondi-Hoyle Accretion - Relativistic Version 

First of all, let us consider the relativistic version of the Bondi-Hoyle ac- 
cretion, i.e., accretion onto a moving nonrotating black hole. The small pa- 
rameter of the problem is again e\ — Woo/coo- In the relativistic case the lin- 
earized stream equation for the stream function $ = <E> [1 — cos 6* + 9)} 
can be written as [Eof 



1 df 
sin 9 89 



£l a 2 N r ?f = 0, (4.57) 
or 



where now 



M 



E 2 



a 2 fi 2 r 2 



(4.58) 



We see that this equation have the same properties as the nonrelativistic 
equation ( |3. 37 ), namely 



• Equation ( 4.57 ) is linear. 

• The values of //, N r , and D should be determined from the unper- 
turbed Schwarzschild metric for a spherically symmetric flow. 

• As ju = fi(r), N r = N r (r), D = D(r), a 2 = 1 - 2M/r, and p = r, one 
can expand the solution in eigenfunctions of the operator Cg. 



But equation (4.57) has another very important property. According to 
(ggg) and Q4.46D, 



D + l 



a 2 



E 2 - u 2 n 2 1 — c 2 



(4.59) 



so the factor a 2 is contained in the all addenda of equation (4.57). Hence, 
equation (4.57) has no singularity at the horizon. In particular, it means 
that it is not necessary to specify any boundary conditions for r = r g . It 
is not surprising because the horizon corresponds to the supersonic region 
which cannot affect the subsonic flow. 

As a result, one can seek the solution of the stream equation in the form 
(for more details see (p^,H) 



$(r,0) = $o[l - cos6» + £i#i(r)sin 2 



(4.60) 
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the equation for the radial function <?i(r) being 



D 



d 2 ffi 
dr 2 



N r 



dffi 
dr 



E 2 - aV 1 - c 2 



— 



0. 



(4.61) 



We see that, as in the nonrelativistic case, the accretion rate is not changed 
in the first order of e±. 

As to the boundary conditions, they are 

1. The regularity condition on the sonic surface. It gives gi{r„) = 0. 

2. At infinity 

9l (r) - K{T) r -, (4.62) 

where (see Table 2) 

1 M 2 r.. 

(4.63) 



2 Ml 



Table 2 



r 


1.01 


1.1 


1.2 


1.333 


1.5 


1.6 


x(r) 


0.49 


0.09 


0.07 


0.044 


0.016 


0.003 


fci(r) 


3.00 


0.56 


0.46 


0.31 


0.12 


0.023 


^in(r) 


-0.74 


-0.09 


-0.03 


0.025 


0.0081 


0.0002 



Equation ( 4.61 ) with boundary conditions 1 and 2 determines the flow 
structure of the Bondi-Hoyle accretion onto a nonrotating black hole. In 
particular, as in the nonrelativistic case, the shape of the sonic surface has 
the form 

r + 1 

2ei— — fci(r) cos< 



i 



(4.64) 



where again fci(r) = g'i(r*)r*. The fci(r) values are given in Table 2. 

But in general, the flow structure remains the same as for the nonrela- 
tivistic accretion (see Fig. 3). Here one can stress only a single interesting 
feature. For r <C r* the radial function gi(r) has the asymptotics 



9l (r) « K in (r) 



-1/2 



(4.65) 



(Ki n is given in Table 2 as well). Hence, for E\ > (M/r*) ' in the vicinity of 
the black hole (i.e., for r < efK? n r^,) the linear approximation is violated, so 



here one should solve the complete nonlinear equation (4.39). Since the sign 
of the coefficient -Km(r) depends on the polytropic index T, the region of 
thickening of the flow lines will either appear on the front side for T > 1.27 
or on the rear side for T < 1.27. However, it takes place in the supersonic 
region which does not affect the subsonic flow for r > . 
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4.2.3 Accretion onto a Slowly Rotating Black Hole 



Let us now consider the accretion of a gas with L = (i.e., a = 2 and 
6 = 2 + 2 — 1 = 3) onto a slowly rotating black hole for which the small 
parameter is 

e 3 = ± « 1. (4.66) 



In this case the metric gik (4.4) differs from the Schwarzschild one by the 
values ~ e\ . One can assume that the thermodynamic functions at infinity, 
Soo and /^oo, remain the same as for spherically symmetric accretion. Hence, 
one can again seek the solution of the stream equation (4.39) in the form 

$(r,0) = $ o [l-cos0 + e£/(r,0)], (4.67) 

where $o is the flux constant corresponding to a spherically symmetrical 
flow. Inserting this form into the stream equation (4.39), we have in the 
first order of e\ 



4) 

where N r is given by (4. 58]). 



(D+ l)sin0 



/' 



d_ 

09 

2 



1 Of 



sin0 89 
M 



eWN r d J- 



E 2 — a 2 fi 2 r 



sin 2 9 cos 9, 



(4.68) 



The properties of equation (4.68) are quite similar to those of equation 



(4.57). In particular, equation (4.68) has no singularity at the horizon. As 
a result, one can show that the flux $(r, 9) can be presented as [^0| 

$(r,0) = $ [(1 - cas6)+e%go{\ - cos<9) + e\g 2 {r) sin 2 9 cos 6}. (4.69) 

Hence, as for the ejection from a slowly rotating star, the flux <&(r, 9) con- 
tains two harmonics, m = and m — 2, the radial function g 2 (r) satisfying 
equation 

M 2 



D 



d^2 
dr 2 



■N, 



d.g 2 



-6 



fj2 



E 2 -oi 2 p? 1 



1 - 2 



M 



The regularity condition on the sonic surface Ng(r* 

^ 1M " 2/ x 

52(n = -x— a ( r *)- 
2 rz 



E 2 — a 2 [i 2 r 

(4.70) 

= gives 

(4.71) 



On the other hand, the condition at infinity (which is the third boundary 
condition after Sqq and gives g 2 (r — > oo) = 0. As a result, the radial 
function g%(r) for r <C r* can be written in the form 



52 M 



M 2 
-G(r)-3- 



(i-3r)/2 



(4.72) 
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where G(r) ~ 1. Hence, at the horizon g2{r g ) ~ (M/r*)( 5_3r )/ 2 , so that 
the disturbance is small (£§<?2( r ) *C 1) everywhere outside a black hole. 
On the other hand, as c^?") < 0, the rotation of a black hole results in the 
concentration of the flow lines in the equatorial plane. Finally, an additional 
consideration demonstrates that p3] 

M 3 

9o = -2— (4.73) 

Thus, the rotation of a black hole diminishes the accretion rate. In reality 
c^q <C 1 and hence M/r* C 1. As a result, the effects of a black hole 
rotation in the vicinity of the sonic surface are actually very small. 

To summarize the last two sections, one can say that the simple zero 
approximation, namely, a spherically symmetric transonic accretion, allows 
us to find analytically 2D structure for very important astrophysical flows. 
A similar approach was used for the construction of an analytical solution 
for 

• accretion of a gas with small angular momentum L (the small param- 



eter e| = (L/Er g ) 2 ) onto a nonrotating black hole |BS , 32|, 



• accretion of a gas with a nonrelativistic temperature (c^ <C 1) and 
without angular momentum (L = 0) onto an arbitrary rotating black 
hole || . 

At present they are the only examples where the analytical solution was 
found. 



4.2.4 Thin Transonic Disk 

As the last example, let us consider the internal 2D structure of a thin tran- 
sonic disk. Such a flow can be realized if the intrinsic angular momentum 
of a gas accreting onto a black hole is large enough (£4 ^> 1) |18|. Here for 
simplicity we consider a nonrotating (Schwarzschild) black hole only (for 
more details see 34 1). 

, B8[, for £4 



According to the standard disk model 35 , 36], [37|, |38| , for £4^1 the 

matter forms a thin balanced disk and performs a nearly circular motion 

1/2 

with keplerian velocity uk(^) ~ (GM/r) . The disk is thin provided that 
the accreting gas temperature is sufficiently low and the disk thickness is 
determined by the pressure gradient 

H^r — . (4.74) 

v K 

Introducing the viscosity parameter ass < 1; relating the stress tensor ££, 
and the pressure as V = ass-P P6L one can obtain 



(4.75) 
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Hence, for c s <C vk the radial velocity v r remains much smaller than both 
the keplerian velocity «k and the velocity of sound c s . 

The General Relativity effects result in two important properties: 

• The absence of stable circular orbits for r < tq = 3r g . 

• The transonic regime of accretion. 

The first point means that the accreting matter passing a marginally stable 
orbit approaches the black hole horizon sufficiently fast, namely, in the 
dynamical time Td ~ [v r (ro)/c]~ 1 ^ 3 r g /c. It is important that such a flow is 
realized in the absence of vis cosity. The second statement results from the 
fact that according to ( |f ,75| ) up to the marginally stable orbit the flow is 
subsonic while at the horizon the flow is to be supersonic. 

Up to now in the majority of works the procedure of vertical averaging 
was used, where the vertical four-velocity us was assumed to be zero p9|. As 



a result, the vertical component of the dynamic force nu b V{,(/iu Q ) in (4.3S) 
was postulated to be inimportant up to horizon. For this reason the disk 
thickness was determined by the pressure gradient even in the supersonic 
region near the black hole |40] |. Here I am going to demonstrate that the 
assumption Ug = is not correct. As in the Bondi accretion, the dynamic 
force is to be important in the vicinity of the sonic surface. 

First of all, let us consider the subsonic region in a close vicinity of the 
marginally stable orbit 7*0 = 3r g where the poloidal velocity u p is much 



smaller than that of sound. Then equation ( 4.39| ) can be significantly sim 
plificd b 
we have 



plified by neglecting the terms proportional to D 1 ~ Up/c^. As a result, 



-M 2 -V k [ — L —V k <S> 



W 1 

a 

64tt 4 / „dE „ AL\ _ , „ds 



> ■ ,„■> \ w E ^~ a L ^ -16tt d nT— = 0. (4.76) 
a 2 w 2 M 2 \ d$ d$ J d<£> 

This equation describing the subsonic flow is elliptical. 

To determine the structure of a two-dimensional subsonic flow (a = 0, 
i.e., 6 = 2 + 3 — = 5) one needs to specify five quantities (three veloc- 
ity components and two thermodynamic functions) on an arbitrary sur- 
face r — ro(8). Naturally, we choose it as the surface of the last sta- 
ble orbit ro = 3r g where ao = a(ro) — y/2/3, w^(ro) = 1/V3, and 
7o = 7( r o) = \/4/3 [l8| . For the sake of simplicity we consider below 
the case where the radial velocity is constant on the surface r = ro and the 
toroidal velocity is exactly equal to u^(ro): 

Uf(r 0) e) = -n , (4.77) 
U e(r ,9) = Buo, (4.78) 
u<p{ro,Q) = 1/V3. (4.79) 



Beskin: 2D Transonic Hydrodynamics 



43 



Here -C \u?\ corresponds to the plane flow at the marginally stable orbit, 
and we introduced the new angular variable 6 = ir/2 — 9 (9disk ~ c ) which 
is counted off from the equator in the vertical direction. Next, we suppose 
that the velocity of sound is also a constant on the surface r = r 

c s (r ,e) = c . (4.80) 



For the polytropic equation of state ( |2.6| ) it means that both the temperature 
Tq = T{tq) and the relativistic enthalpy (1q = /i(fo) are also constant on 



this surface. According to ( 4.75 ), one can find that for a nonrelativistic 



temperature c s <C 1 the small parameter of this problem is 

£ 5 = — ~ a ss co < 1. (4.81) 

CO 

Finally, as the last, fifth boundary condition it is convenient to specify the 
entropy s($). 

Introducing the values Eq = «o7o = and Lo = u l p(ro)ro = V3r g , 

one can write the invariants E(<&) and £($) in the following form 

= fi Q E a = const, (4.82) 
L($) = /x £ocos6 ro . (4.83) 

Here 9 m = 8 m ($) is the angle for which $(r ,6 m ) = $(r, 9). In other 
words, the function S m (r, 9) has the meaning of a theta angle on the last 
stable orbit. This orbit is connected with a given point (r,9) by a line of 
flow $(r, 9) = const. In particular, 9 m (r , 9) = 9. 



First of all, we see that condition E = const (4.82) allows us to rewrite 
equation Q4.76[ ) in a simpler form 

<9 2 $ cos9 d ( 1 _ 2 2 L dL 2 2 2 2 T ds 

9^ + ^59 Vc^9 99 J ~ " ^ d¥ ~ " r C ° S m d¥ " 

(4.84) 



Next, as is shown in Appendix A, for r — ro, the r.h.s. of equation (4.84) 
describes the transverse balance of a pressure gradient and an effective po- 
tential, whereas the l.h.s. corresponds to the dynamic term (vV)v. At 
the marginally stable orbit it is of the order of Uq/cq and may be dropped. 
It is therefore natural to choose the entropy s(<&) from the condition of a 
transverse balance on the surface r = tq 

9 , ds r L dL 

^ C0Se "^ = -«d9V 



where L(9 m ) is determined from the boundary condition (4.83). Thus, we 
have 

s(Q m ) = s(0) - —J ln(cos 9 m ). (4.86) 

6C 
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Owing to ( |4.23| ), one can show that for c s = const relation (4. 
to the standard concentration profile 



i corresponds 



n(ro,6) w noexp 



6c r 



(4.87) 



exercise 



Using relations (2.5)-(2.t), show that the exact expression for n(ro,0) 
is n(ro,6) = n o (cos0) r / 3c o. 



Finally, the definition ( 4.27 ) results in the following connection between 
functions $ and TO 



d<!> = 2na r 2 n(r Q , Q m )u cos6 m d9 r . 



(4.88) 



Hence, due to ( 4.83 ), ( 4.87 ), and ( 4.88 ), the invariant L($) can be directly 



determined fr om th e boundary conditions as well. 

Equation (|4.84|) together with t he bo undary conditions ( 4.83 ), ( 4.86 ), 



and ( 4.88j ) and relationship ( 4.78 ) specifying the derivative d&/dr 
determines the structure of the inviscid subsonic flow inside the marginally 
stable orbit. For example, for a nonrelativistic temperature c s <C 1 we 
obtain 



4 



w 



e 



Here the quantity w, where 



w 2 (r) 



El 



e 2 



a 2 Ll/r 2 



ro 



ro 



(4.89) 



(4.90) 



depending on the radius r only is a poloidal four-velocity of a free particle 
having zero poloidal velocity for r = r$. As we see, w 2 increases very slowly 
when moving away from the last stable orbit. Therefore, the contribution 
of w 2 turns out to be negligibly small and the term may be safely dropped 
in most cases. 

An important conclusion can be drown directly from ( 4.89 ) in which for 
the equatorial plane we have Q r 



O = 0. Assuming u p 



3 = c* and 

neglecting w 2 , we find that the velocity of sound c* on the sonic surface 
r = r* , O = is of the same order of magnitude as the velocity of sound on 
the last stable orbit cq 



r + 4 



Co- 



(4.94) 



Since the entropy s remains constant along the flow lines, the gas concen- 
tration remains approximately constant (n* ~ no) as well. In other words, 
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in agreement with the Bondi accretion, the subsonic flow can be considered 
as incompressible. 

On the other hand, since the density remains almost constant and for 
£5 <§; 1 the radial velocity increases from uq to c* ~ cq, i.e., changes over 
several orders of magnitude, the disk thickness H should change in the same 
proportion owing to the continuity equation (see Fig. 4) 

H{r,) » — H(3r s ). (4.92) 
Co 

As a result, a rapid decrease of the disk thickness should be accompanied 
by the appearance of the vertical component of velocity which also should 



be taken into account in the Euler equation ( 4.38 ). 

Indeed, as one can find analyzing the asymptotics of equation (4.84) j34) , 
in the vicinity of the sonic surface located at 

r* =r - Au 2 Q /3 r , (4.93) 

where the logarithmic factor A = (3/2) 2 ' 3 [ln(co/uo)] ~ 5 — 7, the com- 
ponents of the velocity and the pressure gradient can be presented as 

(4.94) 
(4.95) 
(4.96) 

On the other hand, near the sonic surface the radial scale Sr determining 
the radial derivatives becomes as small as the transverse dimension of a 
disk: Sr w H(r*) w uo^o, so that 

m — (4.97) 
n or 

As a result, both components of the dynamic force 







Uq 


Uf - 








A* 


ul r' 



u 



(4.98) 



_e due ^ c 2 , 6 
r 90 ul r ' 

-4-, (4-99) 
or uf. r 



become of the order of the pressure gradient ( |4.96 ) 



To check our conclusions one can consider the flow structure in the 
vicinity of the sonic surface in more detail. Using again the expansion 
theorem (the smooth transonic flow is analytical at a singular point), one 
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Fig. 4. The structure of the thin accretion disk (actual scale) after passing the 
marginally stable orbit r = 3r g obtained by numerical solving equation (i.84) for 
co = 10 -2 , no = 10 -5 . The solid lines correspond to the range of parameters 
u p/ c ?) < 0-2, where the solution should not differ greatly from the solution of 
the complete equation (4.39). The dashed lines indicate an extrapolation of the 
solution to the sonic-surface region. In the vicinity of the sonic surface the flow 
has a form of an ordinary nozzle 



can write down (cf. ( [2.4g )-(2.4 



n* ( 1 + mh + -T) 3 @' + 



an [ Q + a x hQ + ^a 2 h 2 e + ^oB 3 



(4.100) 
(4.101) 



where h — (r — r*)/r*. Here we assume that all three invariants E, L, and s 
are already given, i.e., i = and 6 = 2 + 0—1 = 1. Hence, as in the planar 
case, the problem needs one extra boundary condition. Now comparing the 
appropriate coefficients in the Bernoulli (4.44) and the full stream equation 
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( 4.43 ), one can obtain neglecting terms ~ Uq/cq 



ao 

ax 
a-2 



9 \ (r+i)/2(r-i) 



M ' 



r + i 



1 - a* 



2.25, 



(i 



(4.102) 

(4.103) 
(4.104) 
(4.105) 

(4.106) 



where a 2 = a 2 (r*) « 2/3. In comparison with the planar case, all the 
coeffi cients are expressed here through the radial logarithmic derivative r\\ 

(HI)- 

Let me stress that it is rather difficult to connect the sonic characteristics 
Vi = ? 7i( r *) with physical boundary conditions on the marginally stable 
orbit r — tq (for this it is necessary to know all the expansion coefficients 
in ( 4.100|) and ( 4.101 )). In particular, it is i mposs ible to form ulate the 
restriction on five boundary conditions ( 4.77 )-( 4.80| ) and ( 4.86| ) resulting 
from the critical condition on the sonic surface. Nevertheless, the estimate 



(4.97) makes us sure that we know the parameter 771 to a high enough 
accuracy. Then, according to (4.102)— ( 4.106) , all the other coefficients can 
be determined exactly. 

The coefficients ( 4.102 )-( 4.106) have clear physical meaning. So, ao 
gives the compression of flow lines: ao = H(tq)/ H(r*). In agreement with 
(4.92) we have ao ~ cq/uo- Further, a\ corresponds to the slope of the flow 
lines with respect to the equatorial plane. As a\ > 0, in a close vicinity of 
the sonic surface the compression of stream lines finishes, so inside the sonic 
radius r < r„ the stream lines diverge. On the other hand, as a\ <C u^ 1 , 
for r = the divergency is still very weak. Hence, in the vicinity of the 
sonic surface the flow has a form of an ordinary nozzle (see Fig. 2). Finally, 
as ai ~ 7/3 ~ bo ~ Uq 2 , one can conclude that the transverse scale of the 
transonic region £f(r#) does the same as the longitudinal one. The latter 
point suggests a very important consequence that the transonic region is 
essentially two-dimensional, and so it is impossible to analyze it within the 
standard one-dimensional approximation. 



Using now the expansions ( 4.100 ) and ( 4.101 ), one can obtain all the 
other physical parameters of the transonic flow. In particular, we have 



I a 2 2 

1 - + -(r - 1) -|e 2 + -(r + 
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2 2 
Ct = C~ 



i + (r - 1) m h + ~(r - 1) ^e 2 - ~(r - i)(r + i)^e 2 



As a result, the shape of the sonic surface u p — c s has the standard parabolic 
form 

h=^ Vl 2 . (4.107) 



Thus, the analysis of the hydrodynamic stream equation (4.39) allows us 
to find a nontrivial structure of the thin transonic disk. As was shown, the 
diminishing disk thickness inevitably leads to an emergence of the vertical 
velocity component of the accreting matter. As a result, the dynamic term 
(vV)v in the vertical balance equation cannot be omitted. In this sense 
the situation is completely analogous to the spherically symmetric Bondi 
accretion for which the contribution of the dynamic term becomes significant 
near the sonic surface and dominant for a supersonic flow. 

However, there is an important difference. For the Bondi accretion the 
dynamic term (vV)v has only one component v r dv r / dr which in the vicinity 
of the sonic surface becomes of the same order of magnitude as both the 
pressure and the gravity gradients. As to the thin accretion disk, both 



components of the dynamic term [(vV)v]g, ( 4.98| ) and ( 4.99| ), become of the 



same order of magnitude as the pressure gradient, the role of the effective 
gravity gradient being unimportant: ~S/eip c s ~ i-e., it is Cq/uq times 

smaller than the leading terms. As a result, the structure of a thin transonic 
disk is quite similar to an ordinary planar nozzle shown in Fig. 2. For 
this reason the critical condition on the sonic surface does not restrict the 
accretion rate. 



5 Conclusion 

Thus, in my lecture I have tried to demonstrate the possibilities and diffi- 
culties of the Grad-Shafranov approach. As we have seen, in some simple 
cases it does allow us to construct the analytical solutions. In particular, 
the Grad-Shafranov approach is very suitable in considering the (analytical) 
properties of the flow in the vicinity of the sonic surface and in determining 
the number of boundary conditions. 

On the other hand, it has been shown that in the general case the regular 
procedure does not exist. The point is that the critical conditions are to be 
specified on singular surfaces which are not known from the very beginning 
and are themselves to be determined from the solution. Moreover, it is 
impossible to extend this approach to nonideal, nonstationary, and non 
axially symmetric flows. For this reason it is not surprising that those 
interested in astrophysics more than mathematics, passed from the Grad- 
Shafranov approach to numerical calculations analyzing absolutely another 
class of equations, namely, the time-dependent ones. The only thing one 
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can wish is not to forget the basic physical results of the Grad-Shafranov 
approach which remain true irrespective on the method of calculation. 

Appendix 



A From Euler to Grad-Shafranov - the Simplest Way 



Let me show how the nonrelativistic version of the stream equation (4.84) 
can be directly derived from the Euler equation. For the sake of simplicity 
we consider here only a case of nonrelativistic velocities and small angles 
= 7r/2 — 9 in the vicinity of the equatorial plane. 
The ^-component of the Euler equation is 

,,2 



dve 
dr 



v e - 



■86 



V r V0 



(A.l) 



Adding v lf> dv v jrd9 to both sides and adding and subtracting v r dv r /rd6 in 
the left-hand side, we get 



dv e 9v r ( v 2 



V r Vg 



r 



dv 



V fl p g . (A.2) 

ro9 m n n 



Using the definition of the Bernoulli integral E — v 2 /2 + w + ip g and the 
thermodynamic relationship dP = m p ndw — nTds, we get for E = const 



dve 
dr 



dv r 
~r~d6 



v r ve 



■ m p ndw - 
, sin 6 d 
rd~9 



K>sin0) + -^-g-. (A.3) 



in. 



rdO' 



It follows from the definition (3.1) that 
1 



<9$ 



1 



r 2Trnr 2 sin 9 d6 ' 9 2nnram9 dr 
Now, assuming n w const (this is the case for a subsonic flow), we obtain 



(A.4) 



1 



47T 2 7 



<9$ 



9 2 $ 
dr 2 



sin 9 d 



1 <9$ 



rVuy sin 



2 d9 \ sm.9 d9 

2 T ds 



sin t 



(A.5) 



Finally, dividing both sides by —{dQ/dff), we get ( 4.84 ). 

Hence, whereas the first term in the l.h.s. of ( 4.84 ) does correspond to 
the component v r dve/dr, and the last term in the r.h.s. (for c s = const) cor- 
responds to the pressure gradient, the role of the term oc LdL/d<& proves to 
be less trivial. It contains both the effective potential gradient and, in fact, 
component vedve/d9. The former is the leading one near the marginally 
stable orbit r w 3r g , whereas the latter becomes important only when ap- 
proaching the sonic surface. 
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